Manifestation of spin-orbit interaction in tunneling between 2D electron layers 
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An influence of spin-orbit interaction on the tunneling between two 2D electron layers is con- 
sidered. Particular attention is addressed to the relation between the contribution of Rashba and 
Dresselhaus types. It is shown that without scattering of the electrons, the tunneling conductance 
can either exhibit resonances at certain voltage values or be substantially suppressed over the whole 
voltage range. The dependence of the conductance on voltage turns out to be very sensitive to 
the relation between Rashba and Dresselhaus contributions even in the absence of magnetic field. 
The elastic scattering broadens the resonances in the first case and restores the conductance to a 
larger magnitude in the latter one. These effects open possibility to determine the parameters of 
spin-orbit interaction and electrons scattering time in tunneling experiments with no necessity of 
external magnetic field. 

PACS numbers: 73.63.Hs, 73.40.Gk, 71.70.Ej 



I. INTRODUCTION 

Spin-orbit interaction (SOI) plays an important role in 
the widely studied spin-related effects and spintronic de- 
vices. In the latter it can be either directly utiUzed to cre- 
ate spatial separation of the spin-polarized charge carries 
or indirectly influence the device performance through 
spin-decoherence time. In 2D structures two kinds of 
SOI are known to be of the most importance, namely 
Rashba and Dresselhaus mechanisms. The first one char- 
acterized by parameter a is due to the structure inver- 
sion asymmetry (SIA) while the second one character- 
ized by (3 is due to the bulk inversion asymmetry (BIA). 
Most brightly both of the contributions reveal themselves 
when the values of a and /3 are comparable. In this case 
a number of interesting effects occur: the electron en- 
ergy spectrum becomes strongly anisotropic^, the elec- 
tron spin relaxation rate becomes dependent on the spin 
orientation in the plane of the quantum well^, a mag- 
netic break-down should be observed in the Shubnilov de 
Haas effect'^. The energy spectra splitting due to SOI 
can be observed in rather well-developed experiments as 
that based on Shubnikov-de Haas effect. However, these 
experiments can hardly tell about the partial contribu- 
tions of the two mechanisms leaving the determination 
of the relation between a and /3 to be a more challenging 
task. At the same time, in some important cases spin 
relaxation time Ts and spin polarization strongly depend 
on the ^ ratio. In this paper we consider the tunneling 
between 2D electron layers, which turns out to be sen- 
sitive to the relation between Rashba and Dresselhaus 
contributions. The specific feature of the tunneling in 
the system under consideration is that the energy and 
in-plane momentum conservation put tight restrictions 
on the tunneling. Without SOI the tunneling conduc- 
tance exhibits delta function-like maximum at zero bias 
broadened by elastic scattering in the layers^, and fluctu- 
ations of the layers width^. Such a behavior was indeed 
observed in a number of experiments^"^. Spin-orbit in- 
teraction splits the electron spectra into two subbands in 



each layer. At that energy and momentum conservation 
can be fulfilled for the tunneling between opposite sub- 
bands of the layers at a finite voltage corresponding to the 
subbands splitting. However, if the parameters of SOI 
are equal for left and right layers, the tunneling remains 
prohibited due to orthogonality of the appropriate spinor 
eigenstates. In^ it was pointed out that this restriction 
can also be eliminated if Rashba parameters are different 
for the two layers. A structure design was proposed^^ 
where exactly opposite values of the Rashba parameters 
result from the built-in electric field in the left layer be- 
ing opposite to that in the right layer. Because the SOI 
of Rashba type is proportional to the electric field, this 



would result in 



-a^, where and are the 



Rashba parameters for the left and right layers respec- 
tively. In this case the peak of the conductance should 
occur at the voltage Uq corresponding to the energy of 
SOI: eUo — ±2akp, where kp is Fermi wavevector. In 
this paper we consider arbitrary Rashba and Dresselhaus 
contributions and show how qualitatively different situa- 
tions can be realized depending on their partial impact. 
In some cases the structure of the electrons eigenstates 
suppresses tunneling at ever voltage. At that the scat- 
tering is important as it restores the features of voltage- 
current characteristic containing information about SOI 
parameters. Finally the parameters a and f3 can be ob- 
tained in the tunneling experiment which unlike other 
spin-related experiments requires neither magnetic field 
nor polarized light. 



II. CALCULATIONS 

We consider two 2D electron layers separated by po- 
tential barrier at zero temperature (see FiglT]). We 
shall consider only one level of size quantization and not 
too narrow barrier so that the electrons wavefunctions 
in the left and right layers overlap weakly. The sys- 
tem can be described by the phenomenological tunneling 
Hamiltonianii^iii: 
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FIG. 1; Energy diagram of two 2D electron layers. 



H = H^ + H(f + Ht, (1) 

where Hq,H^ are the partial Hamiltonians for the left 
and right layers respectively, Ht is the tunneling term. 
With account for the elastic scattering and SOI in the lay- 
ers the partial Hamiltonians and the tunneling term have 
the the following form in representation of secondary 
quantization: 



J2 ^k'^ka^ka ~^ J2 ^kk'^k 
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Here index I is used for the layer designation and can 
take the values / = R for the right layer, Z = L for the 
left layer. By k here and further throughout the paper 
we denote the wavevector aligned parallel to the layers 
planes, a denotes spin polarization and can take the val- 
ues a = ±1/2. e'j, is the energy of an electron in the layer 
I having in-plane wavevector k. It can be expressed as: 



A', 



(3) 



where e — rn being electron's effective mass. Eg and 
A' are the size quantization energy and the energy shift 
due to external voltage for the layer I . We shall also use 
the value A"' defined as A"' = (A'-A'')-h(£o-eo'). Sim- 
ilar notation will be used for spin polarization denoted 
by indices a, a' . The second term in the Hamiltonian 
© ^kk' is matrix element of the scattering operator. 
We consider only elastic scattering. The tunneling term 
Ht in ([2]) is described by the tunneling constant Tkk'aa' , 
which has the meaning of size quantization levels split- 
ting due to the wavefunctions overlap. By lowercase t we 
shall denote the overlap integral itself. Our consideration 
is valid only for the case of weak overlapping, i.e. t <^ 1. 
Parametrically T tsp, where sp is the electrons Fermi 
energy. The term H^q describes the spin-orbit part of 
the Hamiltonian: 

H^so ^a' {(TX k)^ + [3^ {a^k^ - cjyky) , (4) 

where Ci are the Pauli matrices, a', are respectively the 
parameters of Rashba and Dresselhaus interactions for 
the layer I. In the secondary quantization representation: 

HsO = XI '^^V ^ *^^) C^ka^ka'+ih + ik^) C^ka'cVa 
k 

+(3' ~ '''y^ (kx + iky) 4+,cL (5) 
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as^: 



The operator of the tunneling current can be expressed 

(6) 



^ = X! '^kk'aa' {Pkk'aa' 
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We shall assume the case of 



in-plane momentum and the spin projection being con- 
served in the tunneling event so the tunneling constant 

\k'aa' = TSkk'Saa', whcrC S IS the 



,/ has the form Ti 



Cronecker symbol. The tunneling current is then given 
by 



le 
H 



T I dkTr{{pl^)-{pt^)), 



(7) 



where <> denotes the expectation value in quantum- 
mechanical sense. For further calculations it is conve- 
nient to introduce vector operator S'"^,, = I^Oisj = 

{Tr [p'i',,,,,) , Tr (*Tp«;,,,,) }. This vector fully deter- 

mines the current because the latter can be expressed 
through the difference Sqi^ — Sqj^. The time evolution 



of 

^kk' is governed by: 



dS 



kk' 



dt 



(8) 



In the standard way of reasoning we assume adiabatic 
onset of the interaction with characteristic time . We 
will set It; = in the final expression. With this ([8]) turns 
into: 



A(o)n' 

^kk' 



(9) 



Here S'^'^V represents the stationary solution of ^ with- 
out interaction. By interaction here we mean the tunnel- 
ing and the elastic scattering by impurities but not the 
external voltage. The role of the latter is merely shifting 
the layers by ell on the energy scale. From such defined 
interaction it immediately follows that the only non-zero 
elements of -S^^^V' are that with I — I' and k — k' . In 
further abbreviations we will avoid duplication of the in- 
dices i.e. write single I instead of II and k instead of 
kk: 



aim' _ o(0)ir r 
^kk' - *fe ^kk'Ow 

With use of fermion commutation rules 



(10) 



{c^Ck} = {44} = 

{ci4} = 

the calculations performed in a way similar to^^ bring us 



to the following system of equations with respect to S'f. 



II' . 



= (a"' + ihw) Sf + T (Si; - Si) + M(fc)s" 
^, i e' - e - A"' -I- ifiw e - e' - A^'-' + ihw j 



(11) 
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(12) 



where is a known matrix, depending on k and param- 
eters of spin-orbit interaction in the layers. Here we also 
introduced the quadratic forms of the impurities poten- 
tial matrix elements: 



Akk' 
Bkk' = 



- I ^k'k\ 
^k'k^kk' 



(13) 



As ([TT|) and comprise a system of linear integral 
equations these quantities enter the expression ([7|) for 
the current linearly and can be themselves averaged over 
spatial distribution of the impurities. In order to perform 
this averaging we assume the short range potential of 
impurities: 



(14) 



quite naturally, however, it can be similarly traced in the 
technique used in^*^ (see Appendix). For the same rea- 
son the tunneling term should be dropped from (|12p as 
it would give second order in T if substituted into 
(fTT|) . According to (HH) A can be expressed in terms of 
electrons scattering time: 



1 27r /. ^ ,2\ 27r , 



(17) 



where v is the 2D density of states v — 2^j^- By means 
of Fourier transformation on energy variable the system 
PTjl , can be reduced to the system of linear algebraic 
equations. Finally Sif can be expressed as a function of 



mi 



J, . Consequently the current ([7]) becomes a function 

of (j)'^}^^, ^Pfeo'^y -'^'^^ considered case of zero tem- 
perature: 

1 



where 



The averaging immediately shows that the correlators 
{Akk') = A and {Bkk') = B have different parametrical 
dependence on the tunneling transparency t, namely 



(15) 



We emphasize that this result holds for non-correlated 
distribution of the impurities as well as for their strongly 
correlated arrangement such as a thin layer of impurities 
placed in the middle of the barrier. The corresponding 
expressions for these two cases are given below. Index 
'rand' stands for uniform impurities distribution and 'cor' 
for their correlated arrangement in the middle of the bar- 
rier (z = 0): 



B 



rand 



Yin 
W 



dzff{z)f^,{z) r 

dzft (z) 



Yin 
W 



-Qcor 



Yln^ 
W 



f?iO)fl'iO) 



vjn e 

W d 
Ylnl 
W d 



Vins 



W 



w 



d 
1 

W d' 



(16) 



where n and are bulk and surface concentrations of 
the impurities, W is the lateral area of the layers, d is 
the width of the barrier and f{z) is the eigenfunction 
corresponding to the size quantization level, z is coordi- 
nate in the direction normal to the layers planes, z = 
corresponding to the middle of the barrier—. Unlike^'^'' 
and according to p5|) we conclude that the correlator 
{Bkk') has to be neglected as soon as we shall be inter- 
ested in calculating the current within the order of T^. 
In the hereused method of calculation this result appears 



£cr = ± |a' {kx - iky) - /?' {ikx - , 

. Without loss of generality we shall consider the case of 
identical layers and external voltage applied as shown in 
Figffl 
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The calculations can be simplified with account for two 
small parameters: 



C = « 1 

eU ^ 
r]= — < 1 
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(18) 



With (??) calculation yields the following expression for 
the current: 



where 



'^^T'. I I C«) Tr - ded^, 

(19) 
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Parameters a-g are various combinations of the Rashba 
and Dresselhaus parameters of SOI in the layers. Both 
types of SOI are known to be small in real structures so 
that: 



akp <C ep, Pkp ^ £f 



(21) 



This additional assumption together with (|T8|) reduces 
(O to 



2ir 



/ = ^T^^WU J [C^ (sf) + C« (sf)] dip 



(22) 



The integral over (p in can be calculated analytically 
by means of complex variable integration. However, the 
final result for arbitrary a' , /?' is not given here for it is 
rather cumbersome. In the next section some particular 
cases are discussed. 



culations shown below were obtained using the following 
parameters: Fermi energy ep — W meV, spin-orbit split- 
ting was taken to resemble GaAs structures: akp = 0.6 
meV. 



A. No Spin-Orbit Interaction 

In the absence of SOI {a^ = 0, = = 0) the 

energy spectrum for each of the layers forms a paraboloid: 



(23) 



According to our assumptions © , ([7]) , the tunneling takes 
place at: 



(24) 



Both conditions are satisfied only at C/ = so that a 
nonzero external voltage does not produce any current 
despite it produces empty states in one layers aligned 
to the filled states in the other layer (Fig[T]). The mo- 
mentum conservation restriction in (j24|) is weakened if 
the electrons scatter at the impurities. Accordingly, one 
should expect a nonzero tunneling current within a fi- 
nite voltage range in vicinity of zero. For the considered 
case the general formula (|22|) is simplified radically as all 
the parameters (PH)) have zero values. Finally we get the 
well-known resuHi^: 



III. RESULTS AND DISCUSSION 



The obtained general expression p2)) can be simplified 
for a few particular important relations between Rashba 
and Dresselhaus contributions. These calculations reveal 
qualitatively different dependencies of the d.c. tunneling 
current on the applied voltage. The results of the cal- 
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FIG. 2: Tunneling conductance, a: ef 
T = 2 * 10~^^ s; b: same as a, but ak 
as b, but /3 = a; d: same as c, but r = 2 * 10 



10 meV, a = l3 = 0, 
- 0.6 meV; c: same 



2e^T^iyWU- 
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(25) 



The conductance defined as G{U) ~ I/U has Lorentz- 
shaped peak at [/ = turning into delta function at 
T oo. This case is shown in (Figl^la). All the curves in 
Figl2] show the results of the calculations for very weak 
scattering. The corresponding scattering time is taken 
r = 2* 10-"s. 



B. Spin-Orbit Interaction of Rashba type 

The spin-orbit interaction gives qualitatively new op- 
tion for the d.c. conductance to be finite at non-zero volt- 
age. SOI splits the spectra into two subbands. Now an 
electron from the first subband of the left layer can tunnel 
to a state in a second subband of the right layer. Let us 
consider a particular case when only Rashba type of SOI 
interaction exists in the system, its magnitude being the 
same in both layers, i.e. \a^\ — \a^\ = a, (3^ — (3^ — 0. 
In this case the spectra splits into two paraboloid-like 
subbands "inserted" into each other. FigO shows their 
cross-sections for both layers, arrows show spin orienta- 
tion. By applying a certain external voltage Uq 

kp 



2akf 



V2mgf 



the layers can be shifted on the energy scale 



5 



in such a way that the cross-section of the "outer" sub- 
band of the right layer coincides with the "inner" sub- 
band of the left layer (see solid circles in Figl3|). At 
that both conditions are satisfied. However, if the 
spin is taken into account, the interlayer transition can 
still remain forbidden. It happens if the appropriate 
spinor eigenstates involved in the transition are orthogo- 
nal. This very case occurs if = a^, consequently the 
conductance behavior remains the same as that without 
SOI. Contrary, if the Rashba terms are of the opposite 
signs, i.e. — —a^ the spin orientations in the "outer" 
subband of the right layer and the " inner" subband of the 
left layer are the same and the tunneling is allowed at a 
finite voltage but forbidden at [/ = . This situation, 
pointed out in^'^*^ should reveal itself in sharp maxima of 
the conductance at U — ±Uo as shown in Figl^lb. From 
this dependence the value of a can be immediately ex- 
tracted from the position of the peak. Evaluating ([^0]) 
for this case and further the expression we obtain 
the following result for the current: 



/ = 



T 











(26) 



where d = 2akF. The result is in agreement with that 
derived ini^, taken for uncorrelated spatial arrangement 
of the impurities. As we have already noted we do not 
take into account interlayer correlator {Bkk') be- 
cause parametrically it has higher order of tunneling 
overlap integral t than the intralayer correlator {Akk')- 
Therefore the result (|26p is valid for arbitrary degree of 
correlation in spatial distribution of the impurities in the 
system. It is worth noting that the opposite case when 
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FIG. 3: Cross-section of electron energy spectra m the left(a) 
and right (b) layer for the case = -Q!^,/3^ = = 0. 

only Dresselhaus type of SOI exists in the system leads 
to the same results. However, it is rather non-practical 
to study the case of the different Dresselhaus parameters 
in the layers because this type of SOI originates from the 
crystallographic asymmetry and therefore cannot be var- 
ied if the structure composition is fixed. For this case to 
be realized one needs no make the two layers of different 
materials. 



C. Both Rashba and Dresselhaus contributions 

The presence of Dresselhaus term in addition to the 
Rashba interaction can further modify the tunneling con- 
ductance in a non-trivial way. A special case occurs 
if the magnitude of the Dresselhaus term is compara- 
ble to that of the Rashba term. We shall always as- 
sume the Dresselhaus contribution being the same in 
both layers: — (3^ = /3. Let us add the Dressel- 
haus contribution to the previously discussed case so that 
— —a^ = a, a = [3. The corresponding energy spec- 
tra and spin orientations are shown in FigH) Note that 
while the spin orientations in the initial and final states 
are orthogonal for any transition between the layers, the 
spinor eigenstates are not, so that the transitions are al- 
lowed whenever the momentum and energy conservation 
requirement ([M]) is fulfilled. It can be also clearly seen 
from Fig[4]that the condition p4|) . meaning overlap of the 
cross-sections a. and b. occurs only at few points. This 
is unlike the previously discussed case where the over- 
lapping occurred within the whole circular cross-section 
shown by solid lines in FigO One should naturally ex- 
pect the conductance for the case presently discussed to 
be substantially lower. Using we arrive at a rather 
cumbersome expression for the current: 




FIG. 4: Cross-section of electron energy spectra in the left(a) 
and right (b) layer for the case = —a^ = /3. 



where 



G_ (G2_ - G+ {Gl - (52) 



(27) 



G± ^eU±i- 

T 

F± = Gl {Gl - 2S^) . 

Alternatively, for the case of no interaction with impuri- 
ties a precise formula for the transition rate between the 
layers can be obtained by means of Fermi's golden rule. 
We obtained the following expression for the current: 



^ _ 27reT^W ( I ^ Sma'^eU I Sma'^eU ' 



(28) 
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where 



K ^26^ - e^U^ 



Comparing the results obtained from (P7)l and P5|) is an 
additional test for the correctness of (|27p. Both depen- 
dencies are presented in Fig[5] and show a good match. 
The same dependence of conductance on voltage is shown 
in Figl^lc. As can be clearly seen in the figure the con- 
ductance is indeed substantially suppressed in the whole 
voltage range. This is qualitatively different from all pre- 
viously mentioned cases. Furthermore, the role of the 
scattering at impurities appears to be different as well. 
For the considered above cases characterized by reso- 
nance behavior of the conductance, the scattering broad- 
ens the resonances into Lorentz-shape peaks with the 
characteristic width 5 — h/{eT). Contrary, for the last 
case the weakening of momentum conservation, caused 
by the scattering, increases the conductivity and restores 
the manifestation of SOI in its dependence on voltage. 
FiglHd shows this dependence for a shorter scattering 
time r = 2 * lO"""^^. The reason for that is the weakening 
of the momentum conservation requirement due to the 
elastic scattering. One should now consider the overlap 
of the spectra cross-sections the circles in FigU having a 
certain thickness proportional to t~^. This increases the 
number of points at which the overlap occurs and, con- 
sequently, the value of the tunneling current. As the cal- 
culations show, for arbitrary a and /3 the dependence of 
conductance on voltage can exhibit various complicated 
shapes with a number of maxima, being very sensitive to 
the relation between the two contributions. The origin 
of such a sensitivity is the interference of the angular de- 
pendencies of the spinor eigenstates in the layers. A few 
examples of such interference are shown in Fig|Sl a-c. All 
the dependencies shown were calculated for the scatter- 
ing time T — 2-¥ 10~^^ s. FiglHla summarizes the results 
for all previously discussed cases of SOI parameters, i.e. 
no SOI (curve 1), the case ur — —aL,P = (curve 2) 
and an = —ul — P (curve 3). Following the magnitude 
of T all the reasonances are broadenered compared to 
that shown in FiglS] Fig[6lb (curve 2) demonstrates the 
conductiance calculated for the case = ~'^o:r = 
FigElc (curve 2) - for the case az, = \aR = (3. The curve 
1 corresponding to the case of no SOI is also shown in all 
the figures for reference. Despite of a significant scatter- 
ing parameter all the patterns shown in FiglHlremain very 
distinctive. That means that in principle the relation be- 
tween the Rashba and Dresselhaus contributions to SOI 
can be extracted merely from the TV curve measured in 
a proper tunneling experiment. 



IV. SUMMARY 

As we have shown, in the system of two 2D electron lay- 
ers separated by a potential barrier SOI can reveal itself 
in the tunneling current. The difference in spin structure 
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Tunneling conductance calculated for the case = 
/3 and very weak scattering compared to the precise 
result obtained through Fermi's golden rule calculation. 



FIG. 



of eigenstates in the layers results in a sort of interfer- 
ence in the tunneling conductance. The dependence of 
tunneling conductance on voltage appears to be very sen- 
sitive to the parameters of SOI. Thus, we propose a way 
to extract the parameters of SOI and, in particular, the 
relation between Rashba and Dresselhaus contributions 
in the tunneling experiment. We emphasize that unlike 
many other spin-related experiments the manifestation 
of SOI studied in this paper should be observed with- 
out external magnetic field. Our calculations show that 
the interference picture may be well resolved for GaAs 
samples with the scattering times down to ~ 10~^^ s, in 
some special cases the scattering even restores the traces 
of SOI otherwise not seen due to destructive interference. 



ACKNOWLEDGEMENTS 

This work has been supported in part by RFBR, Pres- 
ident of RF support (grant MK-8224.2006.2) and Scien- 
tific Programs of RAS. 



APPENDIX 

In this section we discuss an approach to the calcula- 
tion of the tunneling current based on Green's function 
formalism. As we shall see this approach used in^^ gives 
the same results as that based on the operator motion 
equation ([H]). However, one should be accurate with av- 
eraging over spatial distribution of the impurities. The 
starting point for this calculation is the tunneling Hamil- 
tonian ([T]) in which the coupling between the two layers 
sits merely in the tunneling term Ht (12) • Whenever such 
Hamiltonian is assumed, its part connected to the scat- 
tering at the impurities is the following 
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kk 



„L+ L 
^ka ^k'c 



k,k' 



E 



yR „R+„R 
^kk'^ka Cfe'cri 



(A.l) 



7 



1.0 



0.8 



<i> 

^ 0.6 

CD 

g 0.4 
■a 

8 0.2 



1.0 



0.8 



0) 

o 0.6 

TO 

g 0.4 
§ 0.2 



0.0 



CI . 


— 1 

2 


■ 


3 










3 -2 -1 


1 2 : 


b. 


— 1 


•/\ 


2 - 











7^ Z'), while as was mentioned earlier T tep, hence: 




-2-10 1 
Voltage, mV 



FIG. 6: Tunneling conductance calculated for various param- 
eters of SOI 
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The Hamiltonian (1) can be rewritten as: 
H = Hq + Hy + Hx 



(A.2) 



The calculation of the current can be carried out by 
means of perturbation theory with Hy,HT treated as two 
perturbations. At first step only the tunneling term Ht is 
considered as a perturbation while the rest of the Hamil- 
tonian describes the system treated as unperturbed. The 
tunneling constant T is a small parameter as soon as our 
consideration is restricted to the case of weak tunneling. 
The d.c. tunneling current can be calculated adjusting 
the Kubo formula to the tunneling current The Kubo 
formula treats the perturbation to the leading order so 
that the leading order for the current is T^. We obtain 
the following formula for the d.c. tunneling current ex- 
pressed through the unperturbed Green'a functions: 



Tr / G^V P,e-eU) G^y {p, p, e) dpde. 



(A.3) 

where p is electron lateral momentum, G^y , G^y are the 
components of the unperturbed Green's functions (which 
are matrices in our case). Index V tells that at this 
stage, while only Ht is considered as a perturbation, 
these unperturbed Green's functions include the scatter- 
ing by impurities. At the second step the Green's func- 
tions Gpy are to be expressed in terms of the known 
unperturbed Green's functions Gg of the non-interacting 
2D electron gas. The conventional perturbation theory 
leads to the summation of diagrams of the two tvoes 
□ 



where Vl^, are the matrix elements of the scattering op- 
erator calculated on the electrons eigenfunctions. Note 
that this Hamiltonian does not contain interlayer matrix 
elements V^^, V^^, which would have appeared after 
straightforward secondary quantization of the impurities 
external field p4)) . This is reasonable because these el- 
ements are parametrically small compared to intralayer 
elements V^, and also to the tunneling term Ht- To 
prove the first we recall from (fTB|) that ^ ^/B while 
the intra-layer matrix elements ~ \/A. Taking into 
account ITSl) we have: 



u' 



V' 



t < 1 



The interlayer matrix elements are also parametrically 
small compared to the tunneling term in the Hamilto- 
nian. Indeed, from (HT]) , (dH]) also follows that V" t- 




FIG. 7: Two types of diagrams for calculation of the tunneling 
current 



turbed Green's functions Gg while the crosses correspond 
to the matrix elements V^),, . The ladder diagrams of type 
b. give vertex corrections discussed in^i^. However, in 
our opinion, the ladder diagrams do not give any signif- 
icant impact to the tunneling current. It can be easily 
shown that if the interaction part of the Hamiltonian is 
as given by (jA.ljl . the upper branches of the diagrams 
in Fig[7| (corresponding to Gq) contain only Vj^, while 
the lower branches have only V^^., . Hence the diagrams 
of type a. contain the quadratic forms of the intralayer 
matrix elements Af^^' (I13p while the diagrams of type b. 
contain only interlayer elements Bkk' ■ It follows from 
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(lisp that after averaging over the impurities (with any parameter T and therefore must be omitted in the lead- 
degree of their spatial correlation) the interlayer matrix ing order calculations, 
elements give higher order with regard to the tunneling 
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